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Abstract 

In  this  paper,  we  prove  the  truth  of  the 
following  decomposition  of  the  solutions  to  the 
lossless  wave  equation  in  layered  media:  the 

complete  output  from  a K-layer  media  system, 
which  is  comprised  of  the  superposition  of 
primaries,  secondaries,  tertiaries,  etc.,  can 
be  obtained  from  a single  state  space  model  of 
order  2K  - the  complete  model  - or  from  an 
Infinite  number  of  models,  each  of  order  2K, 
the  output  of  the  first  of  which  is  just  the  pri- 
maries, the  output  of  the  second  of  which  is 
just  the  secondaries,  etc.  This  decomposition 
of  the  solution  to  the  lossless  wave  equation 
into  physically  meaningful  constituents  (i.  e.  , 
primaries,  secondaries,  etc.)  is  called  a 
canonical  Eremmer  Series  decomposition,  after 
Bremmer,  who  in  1951  established  a similar 
dec  ompos  ition. 

In  many  geophysical  situations,  where  reflection 
coefficients  are  quite  small,  the  decomposition 
can  be  truncated  after  secondaries  or  tertiaries; 
hence,  it  also  represents  a way  to  approximate 
the  solution  to  the  wave  equation. 

I.  Introduction 

Nahi  and  Mendel  (Ref.  1)  have  developed  state 
space  equations  for  a model  of  a horizontally 
stratified  nonabsorptive  earth  with  vertically 
travelling  plane  compressional  waves.  Their 
model  is  based  on  a ray  theory  solution  to  the 
lossless  wave  equation  in  layered  media.  As 
such,  it  is  not  limited  just  to  layered  earth 
models;  however,  since  our  interest  is  in  the 
seismic  area,  our  discussions  will  be  in  the 
context  of  such  models.  Mendel  (Ref.  1)  con- 
jectured, and  demonstrated  by  examples,  that 
"the  complete  output  from  a K-layer  media 
system,  which  is  comprised  of  the  superposition 
of  primaries,  secondaries,  tertiaries,  etc.,  can 
be  obtained  from  a single  model  of  order  2K  - 
the  complete  model  - or  from  an  infinite  number 
of  models,  each  of  order  2K  ...,"  the  output  of 
the  first  of  which  is  just  the  primaries,  the 
output  of  the  second  of  which  is  just  the  second- 
aries, etc.  We  refer  to  this  decomposition  of 
the  complete  solution  to  the  lossless  wave 


equation  into  physically  meaningful  components 
(i.  e. , primaries,  secondaries,  etc.)  as  a 
canonical  decomposition. 

In  this  paper,  after  we  review  the  state  space 
model  of  a K-layer  media,  we  prove  the  truth 
of  the  canonical  decomposition  and  Illustrate 
its  use  by  means  of  a simulation  example.  We 
show  that,  whereas  the  complete  model  (i.  e.  , 
the  model  that  generates  primaries  and  all 
multiples  for  the  K-layer  media)  has  both  poles 
and  zeros,  similar  to  an  ARMA  model,  the 
partial  residual  models  (such  as  the  primaries 
model,  which  generates  only  primary  reflections, 
and  the  secondaries  models,  which  generates 
only  the  secondary  reflections)  only  have  zeros, 
similar  to  MA  models.  As  such,  the  partial 
residual  models  are  much  easier  to  compute 
with  than  the  complete  model,  since  they  are 
open-loop  models,  whereas  the  complete  model 
is  a feedback  model.  In  many  geophysical 
situations,  where  reflection  coefficients  are 
quite  small,  the  canonical  decomposition  can  be 
truncated  after  secondaries  or  tertiaries;  hence, 
it  also  represents  a way  to  approximate  the 
complete  solution  to  the  wave  equation. 

After  the  author  had  developed  and  p roved  the 
canonical  decomposition,  he  discovered  that  the 
idea  of  such  a decomposition,  as  well  as  a 
collection  of  integral  equations  for  accomplishing 
it.  were  presented  by  Bremmer  (Ref.  2;  in 
1951.  In  fact,  such  a decomposition  is  referred 
to  as  a Bremmer  Series  (Refs.  3 and  4);  hence, 
we  shall  refer  to  our  canonical  decomposition 
as  a Brecimet  Series  decomposition.  Our 
approach  to  the  development  of  this  decompo- 
sition is  quite  different  from  3remmer's.  We 
shall  have  more  to  say  about  this  in  Section  VII. 

A system  of  K-layered  media  is  depicted  in 
Figure  1.  We  adopt  the  convention  of  calling 
the  layer  below  layer  K the  basement.  Each 
layer  is  characterized  by  its  one  way  travel 
time  , velocity,  V.  , and  normal  Incidence 

reflection  coefficients  r,  (1=1.2 K). 

Additionally,  interface-0  denotes  the  surface 
and  is  characterized  by  reflection  coefficient 
rg  . In  Figure  1,  m(t)  and  y(t)  denote  the  in- 
put (e.  g. , seismic  source  signature  from 


dynamite,  airgun,  etc. ) to  the  layered  earth 
system*  which  is  applied  at  interface-0,  and 
the  output  (i.  e. , ideal  seismogram)  of  the 
system  which  is  observed  at  the  surface,  res- 
pectively. 

As  in  Ref.  5,  we  shall  find  it  convenient  to 
draw  ray  diagrams  with  time  displacement 
along  the  horizontal  axis,  so  that  the  rays 
appear  to  be  at  non-normal  incidence  and  so 
do  not  overlap  one  another. 

Figure  2 depicts  primary  and  multiple  reflec- 
tions for  a 2 layer  media,  and  illustrate  the 
very  complicated  internal  behavior  of  the  two- 
layer  system,  and,  that  y(t)  has  components 
which  occur  at  specific  values  of  t.  The  pri- 
mary reflection  components  of  y(t),  denoted 
y j (t),  are  depicted  in  Figure  3a,  whereas  the 
secondary  reflection  components  of  y(t),  denoted 
y^tt),  are  depicted  in  Figure  3b. 

Observe,  from  this  simple  example,  that  pri- 
mary reflections  result  from  exactly  one  reflec- 
tion off  an  interface  in  a layered  media  system; 
hence,  primary  reflection  components,  y j (t) 
include  only  those  terms  which,  neglecting  the 
effects  of  two-way  transmission  factors,  are 
proportional  to  linear  functions  of  the  reflection 
coefficients  [e.  g.  , rQ(l  - r^)  m(t  - 2Tj)  * rQm(t  - 21^), 
where  * denotes  proportionality].  Observe, 
also,  that  secondary  reflections  result  from 
exactly  three  reflections  off  of  interfaces  within 
a layered  media  system;  hence,  secondary  re- 
flection components,  yjlt),  of  y(t)  include  only 
those  terms  which,  neglecting  the  effects  of 
two-way  transmission  factors,  are  proportional 
to  cubic  functions  of  the  reflection  coefficients. 

In  general,  n-arv  reflections  result  from  ex- 
actly 2n  - 1 reflections  off  of  interfaces  within  a 
layered  media  system;  hence,  n-ary  reflection 
components,  yn(t),  of  y(t)  include  only  those 
terms  which,  neglecting  the  effects  of  two-way 
transmission  factors,  are  proportional  to  pro- 
ducts of  2n  - 1 reflection  coefficients. 

By  thinking  of  y(t)  as  the  superposition  of  its 
constituents  - primaries,  secondaries,  etc.,  - 
we  can  write  it  as 

00 

y(t)  = 2 y.(t)  . (l) 

jel  J 

Equation  (1)  is  the  Bremmer  Series  decompo- 
sition of  y(t).  We  shall  demonstrate  that  such 
a decomposition  is  possible,  and,  shall  provide 


A 

In  a marine  environment,  layer  1 can  be 
taken  as  water.  For  the  purposes  of  the  pre- 
sent paper,  we  assume  that  both  the  source  and 
the  sensor  are  located  at  the  surface.  The 
generalization  of  our  results  to  subsurface 
source  and  sensor  is  straightforward. 


state  space  models  for  obtaining  the  constituents. 

yi(t),y2(t) Our  state  space  models  can 

be  thought  of  as  numerical  procedures  (algo- 
rithms) for  obtaining  the  constituents  of  y(t). 

Let  us  express  each  reflection  coefficient  in 
terms  of  a small  number,  say  t ; that  is  to 
say,  given  r-,  r j, . . . , r^  and  t,  numbers  r^, 
r'j r^  can  be  determined  such  that 

rj  = ‘rj  ’ j = °>  1 k (2) 

From  our  preceding  discussion  on  the  pro- 
portionality of  y^(t)  to  products  of  2i  - 1 reflec- 
tion coefficients,  we  see  that 

y.(t)  = c21'1  y!(t)  . 1*1.2,...  (3) 

which  means  that  y(t)  in  Eq.  (1)  may  be  ex- 
pressed as  the  following  power  series  in  t; 

y(t)  = tyJW  + tV  (t)  + e5y'(t)  + ...  + «Zl"1y!(t)  + . . . 

123  1 (4) 

This  equation  suggests  the  possibility  of  truncating 
the  Bremmer  Series  decomposition  after  a small 
number  of  terms  (unless  the  amplitude  of  yj(t)  is 
comparable  to  (2i-l,  which  is  possible  in  certain 
applications).  It  also  provides  us  with  a method 
for  proving  that  our  partial  residual  state  equation 
models  generate  just  the  primaries,  or  just  the 
secondaries  (these  models  are  described  in  Sections 
IH  and  IV),  etc.  We  shall  prove  that  the  output  of 
our  primaries  model,  yj(t),  is  proportional  to  t, 
that  the  output  of  our  secondaries  model,  y2(t),  is 
proportional  to  t^,  etc. 

In  Section  H we  review  the  state  equation  model 
described  by  Nahi  and  Mendel  (Ref.  1)  for 
generating  the  complete  response.  In  Section 
IH  we  develop  and  prove  the  validity  of  a state 
equation  primaries  model.  In  Section  rv  we 
develop  and  indicate  the  proof  of  validity  of  a 
state  equation  secondaries  model.  A complete 
statement  of  the  canonical  Bremmer  Series  de- 
composition is  given  in  Section  V.  An  example 
which  demonstrates  the  decomposition  for  a 
Bright  Spot  application  is  given  in  Section  VI. 
Relationships  between  our  results  and  Bremmer's, 
as  well  as  conclusions,  are  given  in  Section  VII. 

U.  A State  Equation  Model  for  Complete 

Response 

There  are  (at  least)  two  different  ways  for 
obtaining  a state  equation  model  for  the  com- 
plete response,  y(t),  of  a K-laver  media  system 
(Ref.  1).  The  brief  derivation  of  this  model, 
given  below,  is  based  on  a starting  point  which 
should  be  quite  familar  to  all  exploration  geo- 
physicists. It  is  the  ray  diagram,  depicted  in 
Figure  4,  which  is  taken  from  Ref.  5.  Symbols 
uy  and  d,^  denote  the  upgoing  and  downgoing 
waves  in  the  k-th  layers,  respectively;  and,  we 
adopt  the  convention  that  waves  at  the  top  of  a 
layer  occur  at  present  time,  t.  To  develop  the 
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state  equation  model  we  direct  our  attention  at 
the  intersection  point  of  the  ray  diagram  and 
apply  superposition  to  obtain  the  following 
equations  for  signals  u^(t  - T^)  and  dj  (t), 
which  leave  that  point:  + 

Vt+Tk>  = rvdi(t-Vt'I-V\tl,t|  (5*> 

dk+i(t>  - <5b> 

These  equations  are  applicable  for  k = 1 , 2 

K - 1.  At  the  surface  (Figure  5a),  we  obtain 

y(t)  = rgm(t)  + (1  -rg)Uj(t)  (6) 

dj(t)  = (l  + rQ)m(t)  -rQu1(t)  (7) 


and,  at  the  K4*1  interface,  we  assume  that 
u^+^(t)  = 0,  to  obtain  (Figure  5b) 


K+r 


= 'kV'V 


k'  k' 


K' 


(8) 

(9) 


Signal  y(t)  in  Eq.  (6)  is  the  measurable  system 
output.  Signal  d^+1(t)  is  also  a system  output; 
but,  since  it  cannot  be  measured,  we  shall 
ignore  it  in  following  analyses. 

It  is  convenient  to  group  Eqs.  (5a),  (5b),  (7)  »>nd 
(8)  in  a layer  ordering,  as  follows: 

dl(t)  = -rgUjtt)  + (l  + r0)m(t) 

Vt  + Tl>  * ^dJft-V + (l-r1)u2(t) 


d.(ff  = Cl+*'J.1»«»j.l(*-Tj-1)\ 

- rj-i'Vt) 

u.ft^j)  = rj  dj(t  * Tj, 

+ (1-VVi(t) 


J-2.3. 


, K-l 


dic(t>  * (l  + rK-l,dK-l<t*TK-l)  * rK-luK(t) 


Vt  + V 1 rKdK(t*V 


( 10) 


This  system  of  2K  equations  is  not  in  useful  state 
equation  format,  yet,  since  signals  in  its  left- 
hand  side  occur  at  t and  delayed  times,  and 
signals  on  the  right-hand  side  occur  at  t,  t-Tj.j 
and  t-Tj.  In  order  to  put  Eq.  (10)  into  useful 
state  equation  format,  let 


d.(t)  £ d!(t-T  ) 


(U) 


for  all  j*l,2 K.  Observe,  from  Figure  4. 

that  the  downgoing  states  dj(t)  occur  at  the 
bottom  of  a layer.  Equation  (10)  becomes 

di(t+V  * *rouI(t)  * ir0) m(t) 


Ujlt  + Tj)  = rjd^t)  + (1 -rx)u2(t) 

dj(t + Tj)  = (1+r  >d  (t) 


u.(t  + r.)  = r.d.(t) 

J J J J 

+ (1-rj)Vi(t) 

dK(t  + 'K)  = (1  +rK-l)dK-l(t)  ‘ rK-lUK(t) 
uK(t  + TK)  = rKdK(t) 


By  means  of  transformation  (11)  each  pair  of 
equations  in  (10)  now  only  involves  two  time 
points,  t + Tj  and  t.  Equations  (12)  and  (6) 
together  represent  the  state  equation  model  for 
the  complete  output  y(t);  hence,  they  are  re- 
ferred to  in  the  sequel  as  the  complete  model. 

Comment.  Equations  (5a),  (5b),  (7)  and  (8),  for 
uniform  travel  times  (i.  e.  , Ti=t2=.  . . £ T), 

represent  the  starting  point  for  the  well-known 
transfer  function,  models  of  a K-layer  media 
(Refs.  5,6  and  7,  for  example).  Starting  with 
those  same  equations,  we  have  taken  a path 
different  from  the  one  which  leads  to  transfer 
function  models;  our  path  leads  to  time-domain, 
state  space  equations.  In  the  special  case  when 
T1  =t2=  • • • =tK  = T ’ E9-  (1  2)  can  be  reduced  to  a 
collection  of  finite-difference  equations,  by 
choosing  t = kT.  We  shall  explore  the  impli- 
cations of  taking  this  other  path  on  such  pro- 
blems as  computation  and  deconvolution  in  a 
forthcoming  paper.  ■ 


In  order  to  develop  additional  insight  into  the 
structure  of  our  state  space  model,  it  is  in- 
structive to  portray  that  model  in  signal-flow 
graph  (SFG)  form.  Because  all  important 
features  of  our  model  can  be  observed  from  the 
nature  of  a two-layer  example,  we  shall  pre- 
sent our  SFG  for  that  case,  leaving  generali- 
zations to  the  reader. 


For  purposes  of  constructing  a SFG  it  is  useful 
to  rewrite  Eq.  (12),  for  K = 2,  as  follows: 

dl(t)  = *r0Ul(t- V * (l+r0)m(t- Tj) 


d2(t)  = (l+TjJdjtt-Tj)  - ^ u2(t-T2) 


u2(t)  = r2  d2(t  * t2* 

«l(t|  = ri  di(t- v + (i ' V u2(t ' V 

The  ordering  of  the  state  - dj,d7,uj,Uj  - is  a 
computational  ordering,  because  that  is  pre- 
cisely the  order  in  which  we  would  solve  these 
four  equations.  A SFG  based  on  these  equations 
and  Eq.  (6)  is  depicted  in  Figure  6.  For 
simplicity  in  drawing  the  SFG,  we  do  not  show 
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the  explicit  dependence  of  node  signals  on  time; 
and,  we  use  operators  Zj  and  to  denote  Tj 
and  T£  sec.  time  delays,  respectively  [i.  e. , 
z j m = m(t  - Tj )] . The  feedback  nature  of  the 
layered  media  system  is  clearly  evident  from 
the  SFG.  Three  feedback  loops  are  present; 
the  inner  loop  being  associated  with  layer  2, 
the  other  loops  being  associated  with  layers  1 
and  2. 

The  feedback  nature  of  the  complete  model  in- 
dicates that  it  has  zeros  as  well  as  poles; 
hence,  it  is  similar  to  an  autoregressive  moving 
average  |ARiIA|  model.  It  is  the  feedback 
nature  of  the  complete  model  which  produces 
the  multitude  of  multiple  reflections  (see  Figure 
2,  for  example). 

III.  Primaries  Model 

Conceptually,  the  primary  reflections  (i.  e. , 
primatires)  are  easily  obtained  directly  from 
ray  diagrams  (Figure  3a,  for  example);  however, 
that  approach  does  not  lead  to  a dynamical 
model  which  can  be  associated  with  the  pri- 
maries. In  this  section  we  show  how  to  modify 
the  complete  model  in  a very  simple  manner  so 
as  to  obtain  a state  space  model  which  generates 
only  the  primaries.  We  refer  to  this  model 
either  as  a primary  reflection  model  or  a pri- 
maries model. 

A.  Conjectured  Primaries  Model 

We  conjecture  that  state  equations  for  a pri- 
maries model  are  obtained  from  Eq.  (12)  by 
deleting  the  term  -r.  , u.(t)  in  the  equations  for 
dj(t  + Tj),  j*l,2,...,K.  Our  reason  for  deleting 
these  terms  is  based  on  the  following  observa- 
tion: a multiple  reflection  can  occur  only  when 
an  upgoing  wave  reflects  off  of  the  top  of  a 
layer.  Observe,  in  Figure  4,  that  the  contri- 
bution due  to  “k+j(t)  in  d^+.(t)  represents  just 
such  a reflection;  hence,  we  delete  these  com- 
ponents from  the  complete  model  in  our  pri- 
maries model,  whose  equations  are  given  below. 

dl  l(t+Tl)  * <1  + r0>  m(t) 

Ul,l(t  + V * rl  dl,  l(t)  + (I  ' rl*  “l  2(t) 


dl.J(t*7  3 (1  + rj-l)dl,j-l(t) 
Ul.j(*  + 7 * rjdi,j(t) 


j*2,3, 


+ (1-rj,“l.j  + l(t>  I 
dl,K(t  + V * (1+rK-l)dl,K-l(t) 
ttl.K(t  + V 3 rKdl,K(t) 


y j (*)  * m(t)  + (1  -rQ)  u1(  j(t) 


In  these  equations,  d^j.Uj  :,  and  yj  denote 
primary  downgoing  and  upgoing  state  and  the 
primary  reflection  portion  of  the  complete  out- 
put, respectively. 

At  this  point  Eqs.  (13)  and  (14)  are  conjectural 
in  nature;  i.  e. , that  they  do  give  just  the  pri- 
maries remains  to  be  proven.  We  do  this  in 
Paragraph  3 below.  They  have  been  validated 
by  ray  tracing  arguments  for  a three-layer 
example  in  Ref.  1. 

Before  validating  our  model,  we  wish  to  point 
out  an  important  property  of  the  primaries 
model:  its  SFG  has  no  feedback  paths.  We 
illustrate  this  for  a two-layer  media  in  Figure 
7.  This  property  causes  y j (t)  to  be  comprised 
of  a finite  number  of  terms  - K+l,  for  a K- 
layered  media  - and  enables  y^(t)  to  be  com- 
puted in  a very  straightforward  manner.  The 
open-loop  nature  of  the  primaries  model  in- 
dicates it  only  has  zeros;  hence,  it  is  similar 
to  a moving  average  (MA)  model. 

B.  Validation  of  Primaries  Model 

We  shall  show  that  Vj(t),  as  given  by  Eq.  (14), 
is  proportional  to  eg^(t).  Function  gj(t)  is  net 
equal  to  y ^ (t)  in  Eq.  (3),  because  we  use  pro- 
portionality, rather  than  equality,  arguments  in 
our  proof.  Now,  it  is  possible  that  we  may  not 
have  accounted  for  all  terms  which  are  propor- 
tional to  t in  our  conjectured  primaries  model; 
hence,  we  must  examine  the  residual  between 
y(t)  and  y^(t).  We  shall  show  that  that  residual 
contains  no  terms  which  are  proportional  to  e , 
hence,  y j ( t ) contains  all  the  terms  which  are 
proportional  to  t . and,  therefore,  our  con- 
jectured primaries  model  is  a valid  primaries 
model. 

1.  yj(t)  Proportional  to  c 

We  wish  to  show  that  yj(t),  as  given  by  Eq. 

(14),  is  proportional  to  e g:(t).  To  begin,  we 
rewrite  Eqs.  (13)  and  (14)  in  proportionality 
format,  using  a computational  ordering  of  the 
states,  where  we  assume  (1  ± tr')  f(t)  « f(t).  By 
this  technique,  we  are  absorbing  all  the  two- 
way  transmission  factors,  like  (1  -rf)»(l 
which  appear  in  yj(t).  in  the  proportionality 
arguments.  Equality  arguments  can  also  be 
used;  but.  t)ie  algebra  associated  with  them  is 
more  complicated. 

dl  lftiTl)  * m(t) 


di.j(t+V  * di.j-i(t)  * j = 2<3> 

ttl.K(t  + TK)  * ‘dl,K(t) 


ui.j(t+7  ■ *di.j(t,*ui,j+i(t) 


Yl(t)  “ c m(t)  + u^  j(t)  . (16) 

Obviously,  we  must  show  that  Uj  j(t)  is  linear 
in  c.  From  Eq.  (15),  it  follows  ’that 

dl,2(t  + T2)  * dl,l(t)  * £ »2(‘> 

dl,3(t  + V * d1.2(t)  “ ‘2(t-V  ^ *3(t) 


dl,K(t  + V “ dl,K-l(t)  “ *K-l(t-TK-l)  = *K(t) 


which  demonstrates  that  all  primary  downgoing 
states  are  not  proportional  to  (,  at  all.  From 
Eqs.  (15)  and  (17),  it  also  follows,  that 

“l,K(t  + TK)  “ * dl,K(t)  * * aK(t*TK)  = * bK<t> 

*l.K.l(t  + TK-l)  “ tdl,K-l(t)  + ul.K(t) 

* '•K.l(t-TK.l)MbK(t’TK) 


I t bK_1<t> 


ttl.l(t  + V 'dl,l(t)  + ul,2(t) 

* t m(t  - Tj)  + * b2(t-  r2) 

i t b1(t)  ( 18) 

which  demonstrates  that  all  primary  upgoing 
states  are  proportional  to  s.  Finally,  from  Eq, 
(16)  and  the  last  equation  in  (18),  we  see  that 

y j(t)  * c m(t)  + t bjft)  £ e gjft)  , (19) 

which  was  to  be  shown. 

2.  First  Residual  Model 

We  form  a residual  model  between  our  complete 
model  and  primaries  model.  Let  j(t),  5 ^ # j (t) , 
and  Sj(t)  denote  the  first  residual  upgoing  and 
downgoing  states,  and  first  residual  output. 


respectively 

, where  (j  = 1,2 K) 

Ujft) 

*ui.j(t)  * 

(20a) 

Yj(t)  1 

d.(t) 

" dl.J(t>  • 

( 20b  ) 

and 

*j(t)  * 

y(t) 

- yj(‘)  • 

(21  ) 

Subtract  the  respective  equations  in  Eqs,  (12) 
and  (13),  ar.d  (6)  and  (14)  to  obtain  the  following 
first  residual  state  equations: 


■ij,  + V !<»)  + (1  -fj)  U1>2(t) 

6l.j(t-Tj)  1 (1>rj-l>  Yj-i(t)  \ 


j = 2.3. 
K-l 


Yj(t+7  = rj6i.j(t) 

+ (1-rj)wl,j+l(t)  y 
*l,K(t  + V = <l  + W 4l.K-l(t) 


rK-lUl,K(t)  * rK-lUl,K(t) 
Ul,K(t  + TK>  * rK  *l.K(t) 


»j(t)  = (1 -rQ)  uui(t)  (23) 

By  its  very  definition,  Sj(t)  comprises  all  the 
multiple  reflections  for  a K-layer  media  system. 
Observe  that  the  first  residual  state  equations 
are  driven  by  the  upgoing  states  from  the  pri- 
maries model. 

3.  Sj(t)  Contains  No  Terms  Proportional  to  t 

In  terms  of  proportionality,  Eqs.  (221  and  (23) 
become: 


Si,i(t  + V 


•*ul.l(t)  ■ ,ul.l(t) 


Ul,l(t  + V * ,Sl.l(t)+*l,2(t> 
Yj(t+7*  Yj-i(t,-tui.j(t) 

- euj.lt) 

*l.J(t  + 7 “ * Yj(t,+Ul.j*l(t) 


Y K(t  +’  *k' 


4l.K-l(t)  * •ul.K(t)  * tUl.K(t| 


ll,K(t  + TK>  “ ! 4l.K(t) 


“l.l(t) 


Obviously,  we  must  show  that  Uj  j(t)  contains 
no  terms  which  are  proportional  to  t.  For 
illustrative  purposes,  we  do  this  for  K - 2;  the 
generalization  for  K>2  is  straightforward  and 
can  be  deduced  from  the  topology  of  the  SFC 
for  K>2.  A proportionality  SFC.  in  which 
nodal  signals  are  proportional  to  incoming 
signals  is  depicted  in  Figure  8. 

By  SFG  algebra  or  Mason's  Reduction  Theorem, 
we  obtain  the  following  relationship  between 
Sj(t)  and  Uj_  j(t)  and  Uj  ,(t): 


Yi(t+V  1 *ro-i.i(t)  *roui.i(t> 


s 


I 


f 


I 

i 


L.  | 

i.  . 1 

V i 

' :A 
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ai(t> 


2 2 2 2 2 2 2 
*'  *1*2  ' ' *1(1  + 6 *2 
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t Zj  z2 


1.  1 


(t) 


U1.2(t) 


(26) 


where 


„ , 22  222  22  422 

A = 1 + e z2+e  ZjZ2+t  z^  + e Zjz2 


(27) 


Observe  that  is  common  to  each  term  in  Eq. 
(26).  Additionally,  we  showed  in  Eq.  (18)  that 
Uj  j and  iij  2 are  proportional  to  t;  hence, 
dividing  6 into  the  numerator  operators  in  Eq. 
(26),  we  see  that 
3 


Sj(t) 


s^t)  + other  terms 


(28) 


By  means  of  this  2-layer  example,  we  have 
demonstrated  that  Sj(t)  contains  no  terms  pro- 
portional to  c;  hence,  all  terms  proportional  to 
e have  been  accounted  for  in  yj(t).  This  com- 
pletes the  validation  of  our  primaries  model. 


IV.  Secondaries  Model 


Secondary  reflection  components  of  the  complete 
response  are  much  more  difficult  to  obtain  from 
ray  diagrams  (Figure  3b,  for  example)  than  are 
primary  reflections.  The  main  reason  for  this 
is  that  there  are  many  more  secondaries  than 
primaries.  In  this  section  we  show  how  to 
modify  the  first  residual  model  in  a very  simple 
manner  so  as  to  obtain  a state  space  model 
which  generates  only  the  secondaries.  We  refer 
to  this  model  either  as  a secondary  reflection 
model,  secondaries  model,  or  first  partial 
residual  model. 


A.  Conjectured  Secondaries  Model 

Compare  the  first  residual  3tate  equations  (22) 
with  the  complete  model  state  equations  (12)  to 
observe  that  they  are  structurally  quite  similar. 
Recall  that  in  the  development  of  the  primaries 
model  we  deleted  the  term  -r-  j u.(t)  in  the 
equations  for  d.(t  + r.),  because  that  term  re- 
presented the  S"ignalJ  which  reflects  off  of  an 
interface  and  returns  down  into  the  layered 
media.  We  conjecture,  therefore,  that  the 


term  -r  j _ j ui  j j(t)  in  the  equations  for  4i;(t  + Tj) 
plays  an  analogous  role  to  the  -rj_i  Wj(t)  term; 


hence,  we  delete  the  former  terms  from  the 
first  residual  model  to  obtain  the  following 
secondaries  model  (or.  first  partial  residual 
model); 


d2.1(t*Tl>  " -r0Ul,l(t) 


u2  l(t  + V = rl  d2  l(t)  + (1  "rl>  u2,  2 


(t) 


d2.j(t  + V * (1+rj-l,d2.j-l(t) 


* 'j-l-l.jW 


“z.J^  + V = rjd2.j(t) 


* (1-V  U2,  j+1 


(t> 


J = 2.3, 


, K- 1 


d2.K(t  + TK)  " (1  + rK-l)d2,K-l(t) 


- rK-l  Ul.K(t) 

U2,  K(t  + TK>  = fK  d2,  K(t) 

(29) 

and 

y2(*>  = d-Vu2.i(t) 

(30) 

In  these  equations,  d2  j . u2  j,  and  y2  denote 
secondary  downgoing  and  upgoing  states  and  the 
secondary  reflection  portion  of  the  complete 
output,  respectively. 

Before  discussing  the  validation  of  this  model, 
we  wish  to  point  out  that  the  SFG  of  the 
secondaries  model  has  no  feedback  paths.  W e 
illustrate  this  for  a two-layer  media  in  Figure 
9.  This  property  causes  y2(t)  to  be  comprised 
of  the  requisite  number  of  terms,  and  enables 
y2(t)  to  be  computed  in  a very  straightforward 
manner.  Because  the  secondaries  model  only 
has  zeros,  it  also  is  similar  to  a MA  model. 


B.  Validation  of  Secondaries  Model 


It  can  be  shown  that  y2(t),  as  given  by  Eq.  (30), 
is  proportional  to  e3g2(t).  The  proof  of  thi3 


follows  the  proof  in  Section  III.  B.  1 that  y ^ (t)  is 
proportional  to  e.  First,  show  that  all 


secondary  downgoing  states  are  proportional  to 
e^.  Then  show  that  all  secondary  upgoing  states 
are  proportional  to  s3.  Finally,  from  Eq.  (30), 


show  that  y2(t)“e  g2(t). 


Next,  form  a residual  model  between  our  first 
residual  model  and  secondaries  model.  Let 
J,Z  j(1'  • 42  an<4  s2«)  denote  the  second 

residual  upgoing  and  downgoing  states,  and 
second  residual  output,  respectively,  where 
(j  = l,2 K) 


and 


2.j(t> 

attl.j(t)-U2.j(t> 

(31a) 

2.j(t> 

= 4 j > j(t)  - *2>J(t) 

(31b) 

s,(t) 

= s^t)  - y2(t) 

(32) 

The  second  residual  state  and  output  equations 
can  be  obtained  by  subtracting  the  respective 
equations  in  Eqs.  (22)  and  (29),  and  (23)  and 
(30).  These  equations  will  look  exactly  like  the 
first  residual  state  and  output  equations  in  Eqs. 
(22)  and  (23),  with  4i,j«uitj»  and  »j(t)  re- 
placed by  5 2 j , u2  ■,  and  s-,,  respectively,  and, 
Uj  j (which  forces’  the  first  residual  state 
equations)  replaced  by  u2jj.  The  second  residual 
state  equations  are  driven  by  the  upgoing  states 
from  the  secondaries  model. 


1. 


Final  validation  of  the  secondaries  model  is  to 
show  that  Sjit)  contains  no  terms  which  are 
proportional  to  e2.  The  proof  of  this  follows 
the  proof  in  Section  HI.  B.3  that  Sj(t)  contains 
no  terms  proportional  to  t,  and  is  left  to  the 
reader. 

V.  Canonical  Bremmer  Series  Decomposition 

Our  results  in  the  two  preceding  sections 
suggest  that  a complete  decomposition  of  a 
seismogram  [i.  e. , y (t) ] into  a superposition  of 
outputs  from  primaries,  secondaries,  tertiaries, 
etc.,  models  is  possible.  We  summarize  this 
decomposition  in  the  following: 

Theorem.  The  complete  output,  y(t),  from  a 
K-layer  media  system  can  be  obtained  from  a 
single  model  of  order  2K  - the  complete  model, 
given  by  Eqs.  (12)  and  (6)  - or  from  an  infinite 
number  of  models,  each  of  order  2K,  inter- 
connected as  shown  in  Figure  10.  The  pri- 
maries model  is  defined  by  Eqs.  (13)  and  (14). 
The  n-aries  (or  n-lst  partial  residual)  models 
are  defined  by  the  following: 


dn,l(t  + TI> 


-r0  Un-l,l(t) 


Un,l(t  + V = rldn,l(t)  +(1-rl)un.2(t) 


,(t  + r2)  = (1  +r,  ,)  d 


' j-1'  n,  j-1 

r.  u .(t) 

J-1  n-l,j' 


.(t  + T.)  a r.  d 


j = 2,  3 K-l 


dn.K(t  + V 


+ (1  - r.)  u ,.(t) 
J a,j+l 


= (lfrK-l)dn,K-l(t) 

* rK-l  Vl,K(t) 


,(t  + T„)  = r 


K n.K' 


ya(t)  = (1-ro)un.i(t)  (34) 

for  n=2,  3,...  . Additionally,  y(t)  is  given  by 
Eq.  (1). 

Pr oof:  An  inductive  proof  of  this  theorem  is 

straightforward,  although  details  are  tedious  in 
nature.  We  have  validated  the  primaries  model 
in  Section  HI.  and  have  validated  Eqs.  (33)  and 
(34),  for  n=l,  in  Section  IV.  Assume  the  truth 
of  Eqs.  (33)  and  (34)  for  n = n*  and  show  that 
they  are  valid  for  n=n’*i-l  via  the  inductive  step. 
We  do  not  give  all  the  details,  because  they  are 
so  similar  to  the  ones  in  our  validation  of  the 
secondaries  mcdel  (and  primaries  model).  To 
begin,  ^we  can  show  that  yn*+,(t)  i»  proportional 
to  s2a  +*gn*+.2-  orc*er  *°  <1°  fhis,  we  use  the 
the  (assumed)  fact  that  all  n*-ary  upgoing  states 


[i.  e. , Ujj#  ;(t)]  are  proportional  to  e2a  _1,  to  show 
that  all  (n*+l)-ary  downgoing  states  are  propor- 
tional to  e2n*.  We  then  show  that  all  (n*+l)-ary 
upgoing  states  are  proportional  to  t2n*->-l.  hence, 
/n*+l^  “ t2n*+l  gn*+i(t).  Next>  form  a residual 
model  between  the  n*  residual  model  and  the 
(n*+l)-ary  model.  Let  denote  the  (n*+l) 

residual  output;  i.  e.  , 

V+l(t)  " V‘>  ’ V+l(t)  (35) 

Finally,  show  that  sn*+i(t)  contains  no  terms 
which  are  proportional  to  e2n*+1.  * 

VI.  Example 

To  illustrate  the  canonical  Bremmer  Series  de- 
composition, we  have  chosen  a 3-layer  media 
(Figure  11)  which  can  be  associated  with  a bright 
spot  phenomenon,  because  of  the  thin  low  velocity 
layer  which  is  sandwiched  in  between  the  two  thick 
high  velocity  layers.  We  attribute  no  other  geo- 
logical plausibility  to  this  example.  For  layer  1, 

V i = 7,  500  ft/ sec  and  Sj  = 2.  2 gm/cm2;  for  layer  2. 
V2=5,  500  ft/ sec  and  Oj1!-  6 gm/cm2;  for  layer  3, 
V3=V1  and  p3  = 0j;  for  the  basement,  V4=12,  000 
ft/sec  and  was  approximated  by  the  j-power 
law.  0.  23  V|  (Ref.  8). 

Figure  12  depicts  the  complete  response  as 
well  as  the  primaries,  secondaries  and  some 
of  the  tertiaries  (those  through  2 seconds). 
Occasionally,  a tertiary  and  secondary  will  be 
coincident  so  that  their  sum,  which  appears  in 
the  lower  plot,  looks  different  from  either 
constituent  (see  the  last  spike,  for  example). 
Even  more  interesting  are  the  large  secondaries 
and  tertiaries  which  occur  well  into  the  seismo- 
gram. These  large  components  always  involve 
(at  least  one)  reflection  off  of  the  surface  inter- 
face, where  r^  is  close  to  unity:  however,  that 
in  itself  does  not  explain  such  large  components. 

There  is  a reinforcement  phenomenon  present 
in  a K-layered  media  system,  regardless  of 
the  values  of  the  one-way  travel  times.  Ey 
this  phenomenon,  it  is  possible  for  rays  to 
travel  different  paths  and  yet  reach  the  seismic 
sensor  at  exactly  the  same  time.  In  Ref.  1, 
ray  diagrams  are  given  for  secondary  reflections 
in  a 3-layer  media  system.  There  are  14 
distinct  paths  for  the  secondaries,  but  only  10 
distinct  times  at  which  the  secondaries  will  be 
recorded.  The  four  times  at  which  two 
secondaries  add  together  (i.  e.  , reinforce  one 
another)  are:  t = 4tj  + 2'3  . 4-j  t 2t£  t 2tj  , 

4TJ+4T2  + 2-J,  and  2t  j + 4-2  + 2-Tj  . In  our 
example,  the  reinforced  secondaries  occur  at 
t = 0.  836,  1.  37,  1.406,  and  1.006;  those 
secondaries  are  starred  in  Figure  12.  The 
situation  is  much  the  same  for  the  tertiaries. 

The  large  negative  tertiary  at  t — 1.8  is  due  to 
the  reinforcement  of  a number  of  (as  of  yet 
unknown)  tertiaries.  We  are  presently  studying 


this  reinforcement  phenomenon  to  establish  the 
times  at  which  a reinforcement  occurs  as  well 
as  the  number  of  components  which  add  up  to 
give  the  complete  signal  at  those  times. 

VII.  Conclusions 

We  have  developed  a canonical  Bremmer  Series 
decomposition  of  a complete  seismogram  for  a 
K-layer  media  system.  This  decomposition  is 
depicted  in  Figure  10.  Our  algorithms  for 
computing  the  constituents.  yj(t),  y,(t),  . . . , etc. , 
of  the  complete  seismogram  are  state  space 
primaries,  secondaries etc.  models,  res- 

pectively. These  algorithms  are  summarized 
in  the  theorem  given  in  Section  V.  We  have 
shown  that  each  one  of  our  constituent  models 
only  has  zeros,  whereas,  as  is  well  known,  the 
complete  model  has  both  zeros  and  poles.  This 
property  of  our  constituent  models  makes  them 
easy  to  use  for  computation  purposes;  and,  in 
many  situations  the  Bremmer  Series  can  be 
truncated  after  a relatively  small  number  of 
terms.  Hence,  the  canonical  Bremmer  Series 
decomposition  can  also  be  viewed  as  a com- 
putational technique  for  approximately  solving 
the  lossless  wave  equation. 

We  conclude  by  relating  our  algorithms  for 
computing  the  constituents  yj(t),  y7(t), . . to  those 
given  by  Bremmer  (Ref.  2).  In  Bremmer's 
notation,  Uj  denotes  the  constituent  of  y associated 
with  exactly  j reflections;  hence,  his  Uj.Uj.Ug,  ... 

should  be  associated  with  our  y \ , y 2>  V 3 

Bremmer's  equations  for  his  constituents  are 
given  as  a function  of  spatial  variable  x,  where 
x = 0 at  the  bottom  of  the  layered  media  and  x 
increases  in  the  direction  of  the  basement  to  the 
surface. 

Bremmer's  equations,  from  which  we  can  cal- 
culate u, , Uj,  Ug,  . . . , are  given  by  the  following 
recursive  integrals; 


V*’  = 


k0 

k(x) 


exo-  1 


k(s)  ds 


J2U1(X)  = ' 2[k(x)]ff  a dS 


k (a) 

[k(s)]5  U2i  * 


exp 


{lj*k(j)da} 


( 36a] 


(36b) 


from  knowledge  of  u2 , which  is  calculated 
from  knowledge  of  Uj  ; hence,  Uj  is  calculated 
from  knowledge  of  uj  ; Ug  is  calculated  from 
knowledge  of  u3  , etc.  Once  the  primaries,  Uj, 
have  been  computed,  we  can  combine  Eqs.  (36b) 
and  (36c)  so  that  we  can  view  u2i+1  as  being 
calculated  from  knowledge  of  u2^_i  (i=l,2,...). 
This  dichotomy  between  the  algorithm  used  for 
computing  the  primaries  and  the  algorithm  used 
for  computing  all  other  constituents  also  appears 
in  our  Section  V algorithm. 

In  our  development,  we  have  only  been  interested 
in  the  surface  seismogram.  y(t).  In  Bremmer's 
development,  the  surface  is  at  x = 0;  hence,  his 
constituents  Uj(x),  Uj(x),  . . . , are  related  to  our 
constituents  yjftl.yjft), . . . , by  the  relationship 

that  y;(t)  = U2i_i  (*=0).  j*l,2 Just  as  he 

has  related  u2i+1  to  u2^_j  , we  have  related  our 
upgoing  state  vectors  (see  Figure  10)  iin  to 
un_j  , so  that  yn(t)  can  be  computed  using 
knowledge  of  un_j  . His  recursive  relationships 
involve  the  constituents  directly,  whereas  our 
recursive  relationships  involve  the  states  for 
successive  partial  residual  models. 

It  would  be  interesting  to  explore  the  relation- 
ships between  Bremmer's  algorithm  and  ours 
in  more  detail.  This  remains  to  be  done. 
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where  i * 0,  1,  2,  . . 


(36c) 

denotes  the  constant  wave 


number  of  the  homogeneous  half  space  x<0;  and 
k(x)  represents  the  variable  wave  number  in  an 


These  equations  were  derived  by  replacing 
the  inhomogeneous  space  x>0  by  a set  of  homo- 
geneous layers  0 <x  <Xj  , Xj  <x  <x2  , Xi  <x  <x2  , . . 
with  the  successive  constant  wave  numbers 

kj.k-.kj At  a later  point,  Bremmer 

passes  from  this  discontinuous  medium  to  a 
continuously  changing  one  by  making  the  thick- 


inhomogeneous  space  x>0.  Observe  that  Uj  is  ness  Ax3  = xs-x3_|  of  the  various  layers  in- 

calculated  from  knowledge  of  Uq;  Uj  is  calculated  finitely  small. 
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Figure  1 System  of  K-layered  media. 


Figure  2 Two  layer  example. 


Figure  3 Primaries  (a)  and  secondaries  (b)  for 
a two  layer  media.  In  (b),  dashed 
lines  denote  the  primaries  which 
cause  respective  components  of  the 
secondaries. 
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Figure  4 


Reflected  and  transmitted  waves  at 
interface  k.  From  Equation  (11), 


v*- V = • 


Figure  5 Reflected  and  transmitted  waves  at 
(a)  surface  and  (b)  interface  K. 
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Figure  9 SFG  for  secondaries  model,  for 
K = 2. 
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Figure  6 SFG  based  on  computational  ordering 
of  nodes  (states)  for  K = 2. 


Figure  7 SFG  for  primaries  model,  for  K - c. 


Figure  3 Proportionality  SFG  for  first  residual 
model,  for  K = 2.  This  SFG  is 
obtained  using  Equations  (24)  and  (25) 
for  K = 2. 


Figure  10  Canonical  Bremmer  Series  decom- 
position of  a seismogram  signal, 
y(t).  Vector  iin  denotes  the 
collection  of  K upgoing  states  for 
the  n-aries  model. 
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Figure  11  Three  layer  bright  spot  example. 


